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The Josephson ring modulator (JRM) is a device, based on Josephson tunnel junctions, capable of 
performing non-degenerate mixing in the microwave regime without losses. The generic scattering 
matrix of the device is calculated by solving coupled quantum Langevin equations. Its form shows 
that the device can achieve quantum-limited noise performance both as an amplifier and a mixer. 
Fundamental limitations on simultaneous optimization of performance metrics like gain, bandwidth 
and dynamic range (including the effect of pump depletion) are discussed. We also present three 
possible integrations of the JRM as the active medium in a different electromagnetic environment. 
The resulting circuits, named Josephson parametric converters (JPC), are discussed in detail, and 
experimental data on their dynamic range are found to be in good agreement with theoretical 
predictions. We also discuss future prospects and requisite optimization of JPC as a preamplifier 
for qubit readout applications. 
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The photon energy of microwave radiation in the band 
from 4 — 8 GHz (^8 — 4 cm wavelength) is approxi- 
mately 10 5 smaller than that of the visible light. Yet, 
at a temperature 10 4 smaller than room temperature, 
now routinely achievable with a dilution refrigerator, it 
is now possible to resolve the energy of single microwave 
photons 1 . There are three advantages of single photon 
microwave electronics when compared with quantum op- 
tics. First, signal shapes at carrier frequencies of a few 
GHz with a relative bandwidth of few percent can be 
controlled with much greater relative precision than their 
equivalent at a few hundreds of THz. This is partly due to 
the fact that microwave generators have more short term 
stability than lasers, but also because microwave com- 
ponents are mechanically very stable, particularly when 
cooled, compared with traditional optical components. 
Second, in single photon microwave electronics, the on- 
chip circuitry can be well in the lumped element regime, 
and spatial mode structure can be controlled more thor- 
oughly and more reliably than in the optical domain. Fi- 
nally, there exists a simple, robust non-dissipative com- 
ponent, the Josephson tunnel junction (JJ), whose non- 
linearity can be ultra-strong even at the single photon 
leveP. Many quantum signal processing functions have 
been realized using JJs, both digital and analog, and this 
short review will not attempt to describe all of them. We 
will focus on analog Josephson devices pumped with a 
microwave tone. They recently led to micr owave ampli- 
fiers working at the single photon leveP^J. These novel 
devices have taken the work pioneered by B. Yurke at 
Bell labs 25 years agcEHU to the point where actual ex- 
periments can be performed using Josephson amplifiers 
as the first link in the chain of measurement^" 1 ^. 

In this paper, we address one particular subclass of 
analog signal processing devices based on Josephson tun- 
nel junction, namely those performing non-degenerate 
three-wave mixing. Examples are Josephson circuits 
based on the Josephson ring modulator * 11 * 12 * which we 
will describe below. The Hamiltonian of such a device is 



of the form 
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+~ {JC X X 2 + JC Y Y 2 + K, Z Z 2 ) + KXYZ, (1) 

where (X, Y, Z) and (Px, Py, Pz) are the general- 
ized position and momentum variables for the three in- 
dependent oscillators, A4x,y,z and ICx,y,z represent the 
"mass" and "spring constant" of the relevant oscillator 
(see table I), and K is the three- wave mixing constant 
which governs the non-linearity of the system. We will 
discuss later how such simple minimal non-linear term 
can arise. The classical equation of motions for the stand- 
ing waves in such a device are symmetric and are given 
by: 

X + 7 a Y + lo 2 X + K'YZ = x(t) cos Lo a t, (2) 
Y + lb Y + ulY + K'XZ = y{t) cos oj b t, (3) 
Z + j c Z + lj 2 Z + K'XY = z(t)cosoj c t, (4) 

where K' = K/A4 (we assume, for simplicity, equal 

masses Mx,y,z = M) and uj a ,b,c = \/K-x,y,z / 'M. are 
the angular resonant frequencies of the three coordinates 
satisfying 

LU a < LO b <CJ c =L> a + UJ b . (5) 

We also suppose the oscillators are well in the under- 
damped regime 

7a < w a , (6) 
lb < w&) (7) 
7c < w c , (8) 

a sufficient but not strictly necessary hypothesis, which 
has the principal merit of keeping the problem analyti- 
cally soluble under the conditions of interest. It is worth 
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X (position) 


P (momentum) 


A4 (mass) 


K, (spring constant) 


<E> (flux) 


Q 


C (capacitance) 


L" 1 


Q (charge) 




L (inductance) 


C" 1 



TABLE I: Generalized variables and parameters for the system of oscillators described by Eq. (1). The variables and parameters 
listed in the first line apply to the case of a mechanical oscillator, whereas the ones listed in the second and third lines are 
adapted for describing an LC oscillator with parallel and series dissipations, respectively, as shown in Fig. 1 and Fig. 16 (top 
panel) for the parallel case, and Fig. 16 (bottom panel) for the series case. 



noting that the system is non-degenerate both spatially 
and temporally. On the other hand, it is important to 
suppose that the envelope functions x(t), y (t) and z (t) 
of the drive signals are supposed to be slow compared to 
the respective drive frequencies w — u> a ^> j a + 75. 

The equations (2|4) must be contrasted with that of a 
degenerate three- wave mixing device for which two cases 
are possible. In the first case, where the Y and Z de- 
grees of freedom have merged into a single oscillator, the 
Hamiltonian has a non-linear term of the form KXZ 2 
and the equations read: 

X + j a X + u? a X + K'Z 2 = x (t) cos ojJ, (9) 
Z + j c Z + u*Z + 2K' ZX = z(t)cosuj c t. (10) 

This is the case of electromechanical resonators^ in 
which one of the capacitance plates of a microwave oscil- 
lator (Z) is itself the mass of a mechanical resonator (X). 
There uj c 3> uj a , and pumping the microwave oscillator in 
the vicinity of oj c — to a leads to cooling of the mechanical 
oscillator provided 7 C 3> 7a- In the second case, it is the 
X and the Y degrees of freedom that merge into a single 
oscillator, leading to a non-linear term in the Hamilto- 
nian of the form KX 2 Z. The equations then read 

X + -y a X + u)lX + 2K'XZ = x(t)cosui a t, (11) 
Z + lc Z + uj 2 Z + K'X 2 = z (t) cos uj c t (12) 
and we have now 
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(13) 



This case is implemented in Josephson circuits as a dc- 
SQUID whose flux is driven by a microwave oscillating 
signal at twice the plasma frequency of the SQUIE 14 . 
When z (t) — z d > K'X 2 (so-called "stiff' or "non- 
depleted" pump condition), the system of equations 
( 11|12 1 reduces to the parametrically driven oscillator 
equation 



X + la X 



K'z d 

7cW c 



sin (u! c t) 



X = x (t) cos UJ a t. 



(14) 

Note that there is, in addition to the parametric drive 
on the left hand side, a small perturbing drive signal 
x (t) cosui a t on the right hand side. The theory of the 
degenerate parametric amplifier starts with this latter 



equation, the term K * d sin (oj c t) corresponding to the 
pump and x it) cos cu a t corresponding to the input signal. 
The output signal is obtained from a combination of the 
loss term ^ a X and the input signal. 

In the context of Josephson devices, another route to 
the effective parametric oscillator of equation ( 14 1 can be 
obtained by a driven, Duffing-type oscillatoi f la | X0 1 . This 
system (Josephson bifurcation amplifier) has only one 
spatial mode and quartic non-linearity, 

X + la X + u? a X - XX 3 = [z d + x (*)] cosLo d t. (15) 

Driven by a strong tone z d cos ui d t in the vicinity of the 
bifurcation occurring at 
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u d = w a — 7 a , 



z d 



(16) 
(17) 



it will lead to an equation of the form ( 14 ) for small devia- 



tions around the steady-state solution. It will, therefore, 
amplify the small drive modulation signal x (t) of equa- 
tion (151 Q~7]. Similar amplifying effects can be found in 
pumped superconducting microwave resonators without 
Josephson junction 3 18 " 20 [ 

In the following section, we will treat Eqs. ([2]|4]) 
using input-output theorySB and obtain the quantum- 
mechanical scattering matrix of the signal and idler am- 
plitudes in the stiff-pump approximation. This allows us 
to find the photon gain of the device in its photon am- 
plifier mode as a function of the pump amplitude, and 
the corresponding reduction of bandwidth. We then dis- 
cuss the implementation of the device using a ring of 
four Josephson junctions flux-biased at half-quantum in 
Sec. II. It is the non-dissipative analogue of the semi- 
conductor diode ring modulator^. In Sec. Ill, we treat 
the finite amplitude of signals and establish useful rela- 
tions between the dynamic range, gain and bandwidth. 
In Sec. IV we introduce the Josephson parametric con- 
verter (JPC) as an example of a non-degenerate, three- 
wave mixing device operating at the quantum limit. We 
present three different realizations schemes for the JPC 
and point out their practical advantages and limitations. 
In Sec. V we present experimental results for different 
JPC devices and compare the data with the maximum 
bounds predicted by theory. We follow this with a discus- 
sion, in Sec. VI, of general requirements for an amplifier 
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to meet the needs of qubit readout and how the maximum 
input power of the device can be increased by two orders 
of magnitude beyond typical values achieved nowadays. 
We conclude with a brief summary of our results in Sec. 
VII. 



I. INPUT-OUTPUT TREATMENT OF A 
GENERIC NON-DEGENERATE, THREE- WAVE 
MIXING DEVICE 



The three oscillators of Eqs. (|2|4j) correspond to three 
quantum LC oscillators coupled by a non-linear, trilin- 
ear mutual inductance, whose mechanism we will discuss 
in the next section. They are fed by transmission lines 
which carry excitations both into and out of the oscilla- 
tors, as shown on Fig. [I] The Hamiltonian of the system 
is (leaving out the transmission lines for the moment), 




FIG. 1: General non-degenerate three-wave mixing device 
consisting of three LC oscillators coupled by a non-linear 
medium, giving a trilinear term in the Hamiltonian of the 
form isT$ a $i,$ c where the fluxes $ a ,t,c are those of the in- 
ductors. Each oscillator is fed by a transmission line with 
characteristic impedance R a ,b,c- 



— - = uj a a) a + Uffb^b + w c c^c 
n 

+ g 3 (a + a^) (b + tf) (c+c*) , 



(18) 



where a, b and c are the annihilation operators associ- 
ated with each of the three degrees of freedom. Their 
associated angular frequencies are given in terms of the 
inductances and capacitances as 
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yj L a ,b,cC a ,b,i 



(19) 



The bosonic operators of different modes (a, b, c) com- 
mute with each other and those associated with the same 
mode satisfy the usual commutation relations of the form 



[a,a f ] = 1. 



(20) 



The link between the mode amplitude such as X, which 
represents the flux through the inductance of the oscil- 
lator, and a quantum operator such as a can be written 

as, 

X = X ZPF (a + a^), (21) 
where "ZPF" stands for "zero-point fluctuations" and 



X 



ZPF 



z n = 




(22) 
(23) 



the last equation defining the impedance of the oscilla- 
tor, equal to the modulus of the impedance on resonance 
of either the inductance or the capacitance. The link 
between K and 33 is therefore 



Hg 3 = KX 



ZPFyZPF ^ZPF 



(24) 



We now work in the framework of Rotating Wave Ap- 
proximation (RWA), in which we only keep terms com- 
muting with the total photon number 



— aJaO? a+Ubb^ b+uicC 1 c+g3 (a^b^c + abc*) . (25) 



Treating in RWA the coupling of each oscillator with a 
transmission line carrying waves in and out of the oscil- 
lator (see Appendix for complements of the next 6 equa- 
tions), one arrives at three coupled quantum Langevin 
equations for a(t), b (t) and c (t) : 



d . . , t 7a 

— a = —iio n a — iq-iw c a 

At 2 

d u ■ h i „ t T fc h 
— b = —tLUhb — zo^a'c - 

df 2 
d 



/7a~« m ' 



df 



mab - —c 



(*), 
Tbb™ (t) , 
Tcc in (t) , 



(26) 



In these equations, the second term in the right hand 
side corresponds to the non-linear term producing photon 
conversion. The third term says that photons introduced 
in one resonator leave with a rate 



Ta,o,< 



Ra,b.< 



(27) 



with the resistances R a ,b,c denoting the characteristic 
impedances of the transmission lines. Finally, in the 
fourth term of the Langevin equations, the input fields 
such as a m (t) correspond to the negative frequency com- 
ponent of the drive terms in the classical equations. They 
obey the relation 



a in (i) = 



1 



+ OO 



(28) 







where a m [tj] are the usual field operators obeying the 
commutation relations 



[a in H 



sgn 



S(oj + uj') (29) 
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in which u denotes a frequency that can be either positive 
or negative. The transmission lines thus both damp and 
drive the oscillators. The incoming field operator treats 
the drive signals and the Nyquist equilibrium noise of the 
reservoir on the same footing. Photon spectral densities 
A/" ln [w] of the incoming fields, introduced by relations of 
the form 



({a i »[u;],a i »[u;'}}) = 2Af™ 



6(u> + u>'), (30) 



have the value 

M 



sgnH coth 



huj 



2 \2k B T / 

+ 2irP™ [6 (w - Ul ) + S(u + wi)] 



(31) 



where P™ is the photon flux of the incoming drive sig- 
nal at angular frequency u)\ (in units of photons per unit 
time) and T is the temperature of the electromagnetic ex- 
citations of the line. Note that the dimensionless function 
Af" 1 [w] is defined for both positive and negative frequen- 
cies. It is symmetric Af™ [uj] = M™ [— uj] and its value at 
frequency \lo\ represents the average number of photons 
per unit time per unit bandwidth in the incoming signal, 
which in the high temperature limit is fc^T/ (fi.|w|). It 
includes the \ contribution of zero-point quantum noise. 

It is worth insisting that we treat the non-linear cou- 
pling strength as a perturbation compared with the in- 
fluence of the reservoirs, treated themselves as a pertur- 
bation compared with the Hamiltonian of the oscillators: 



ff3 < 7a, lb < 7c < W OJ LO b < u c 



(32) 



In general, only one strong drive tone is applied to one of 
the resonators and is called the "pump" . Two cases must 
then be distinguished at this stage, as shown in Fig. [2j 

Case 1 (amplification and frequency conversion with 
photon gain): the pump tone is applied to the c res- 
onator. The device is usually used as an amplifier^li^. It 
can also be used as a two-mode squeezer^. 

Case 2 (noiseless frequency conversion without pho- 
ton gain): the pump tone is applied to either the a or 
b resonatoi^. The device is useful as a noiseless up- 
and down-converter and can perform dynamical cooling 
of the lowest energy oscillator, transferring its spurious 
excitations to the highest frequency one, which is more 
easily void of any excitations and plays the role of a cold 



CASE 1 : 
AMPLIFICATION 
AND CONVERSION 
WITH PHOTON GAIN 
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FIG. 2: Characteristic frequency landscape of non-degenerate 
three-wave mixing devices. Three separate oscillators have 
resonant frequencies ui a < u)b < oj c — u) a + uib ■ They are fed 
by transmission lines, giving them a full linewidth at half- 
maximum 7 a , 7b and 7 C respectively. The non-linear coupling 
strength, expressed in photon amplitude language, is much 
smaller than these linewidths. The device can be pumped 
at uj c and operates then as a phase-preserving amplifier with 
photon gain for frequencies uj a and ujb (top), or it can be 
pumped at one of the two lower frequencies uj a or uib and 
operates then as a noiseless frequency converter or dynamical 
cooler, upconverting signals into oscillator at uj c (bottom). In 
this figure, the spectral density of weak signal corresponds 
to thin arrows whereas the spectral densities of pump signals 
corresponds to thick arrows. 



the signal and idler photons at uj a and ujf,. For solving 
the quantum Langevin equations, we replace the pumped 
oscillator annihilation operator c by its average value in 
the coherent state produced by the pump as 



c(t)^(c(t)) = V^e- i ^*+^. 



(35) 



A. Photon gain (case 1) 



The Langevin equations can then be transformed into the 
linear equations (see equation (1781 of Appendix) 



We will first suppose that the pump is "stiff" , namely 

|(5 in )| 2 »l (33) 

7c>7a,7f> (34) 

This means that the pump tone will not be easily de- 
pleted despite the fact that its photons are converted into 



oi 



or 



C&out 



O- W a b e 



a p — iuj c t 



a'" 
feint 



(36) 
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where 



d_ 

7o,b 

2 : 



r a ,b 



a,b /=— _i<A r^a,b 



(37) 
(38) 

(39) 



After a Fourier transform, we obtain in the frequency 
domain, a simpler relation 



/i a [uii] +ig% 
-Wl* K [wa] 



a out 

6 out [-W2] 



-w 2 ] 



where 



(40) 



(41) 



and the signal and idler angular frequencies cj% and io 2 
are both positive, satisfying the relationship 



LUl + LJ 2 = 0J C . 



(42) 



The scattering matrix of the device for small signals is 
defined by 



(43) 



a out 




T aa Sab 




a in [+ui] 


6 out [-W2] 




Sba r b b 




_ b in [-u> 2 ] _ 



It can be computed from Eq. ( 40 ) and one finds 



Xa^Xb 1 * + W 



(44) 

(45) 
(46) 
(47) 



where the x's are the bare response functions of modes 
a and b (whose inverses depend linearly on the signal 
frequency) 



' aa ~ 


-1 -1* 

Xa X b - 


- \P\ 2 


rbb = 


Xa^Xb 1 + 


\P\ 2 


-1 -1* 
Xa X h - 


\ P r 


Sab = 


-2ip 




-1 -1* 

Xa Xb - 


\ P r 


Sba = 


2ip* 




-1 -1* 

Xa X b - 


\ P r 



Xa 



l-i 



LOl - LO a 



X b = l - l —n ' 

1 b 

and p is the dimensionless pump amplitude 

-i<f> 

'i. ; \/ ' ' < 

P 



VTaTb 



(48) 
(49) 

(50) 



Note that the matrix in Eq. ( 43 ) has unity determinant 
and the property 



nb 



I Sob I = !> 
2 



Sba 



1. 



(51) 
(52) 



For zero frequency detuning, i.e. x a 1 = Xb 1 = 1) the 
scattering matrix displays a very simple form 



cosh To 
+ie +i<t ' sinhr 



-ie ^sinhro 
cosh t 



(53) 



where tanh(ro/2) = \p\. The zero frequency detuning 
power gain Go is given by 



G = (cosh r ) 



l-\pf 



(54) 



For non-zero detuning, the scattering matrix acquires 
extra phase factors but the minimal scattering matrix 
for a quantum-limited phase-preserving amplifier repre- 
sented in Fig. [3] still describes the device. The gain Go 



port b 




FIG. 3: An amplifier reaching the quantum limit must have 
a minimal scattering matrix, with the signal in port a being 
reflected with amplitude gain G 1 ^ 2 while the signal in port b 
is phase-conjugated and transmitted to port a with amplitude 
gain (G — 1) 1//2 . This can be realized in case 1 of Fig. [2] 



diverges as \p\ — > 1 , i.e. when the photon number n c in 
the pump resonator reaches the critical number given by 



r Q r b 
M 2 



(55) 



a result that is common to all forms of parametric am- 
plification. Increasing the pump power beyond the criti- 
cal power yielding leads to the parametric oscillation 
regime. This phenomenon is beyond the scope of our 
simple analysis and cannot be described by our starting 
equations, since higher order non-linearities of the sys- 
tem need to be precisely modelled if the saturation of 
the oscillation is to be accounted for. 
Introducing the detuning 



Aw = LOl 



■ UJ b - L0 2 , 



(56) 



(> 



we can give a useful expression for the gain as a function 
of frequency as 



G(Aw) 



Go 



1 



7G " 



(57) 



which shows that in the limit of large gain, the response 
of the amplifier for both the signal and idler port is 
Lorentzian with a bandwidth given by 



B = 2 7 G, 



-1/2 _ 27 a 7bG 



1/2 



la + 7b 



(58) 



The product of the maximal amplitude gain times the 
bandwidth is thus constant and is given by the harmonic 
average of the oscillator bandwidths. Another interest- 
ing prediction of the scattering matrix is the two-mode 
squeezing function of the device demonstrated in Ref. [S5J 



B. Conversion without photon gain (case 2) 

The case of conversion without photon gain can be 
treated along the same line as in the previous subsection, 
where scattering takes place between c and a or c and 
b modes. Without loss of generality we assume that the 
pump is applied to the intermediate frequency resonance. 
In this case the scattering matrix reads 



scattering matrix is now unitary (conservation of total 
number of photons) and satisfies the following relations: 



\r aa \ 2 + \t ac \ 2 = 1, 
\r cc \ 2 + \t ca \ 2 = 1. 



(63) 
(64) 



For zero frequency detuning, i.e. \ a 1 
scattering matrix can be written as 



cos ro e l ^smTo 
cos To 



sin r 



(65) 



which corresponds to replacing the parameter To by ir^ 
or \p\ by i \p\ in the scattering matrix (53). A scattering 



representation of the two-port device in conversion mode 
is shown in Fig. |4j In this mode the device operates as a 
beam splitter, the only difference being that the photons 
in different arms have different frequencies®!. Full con- 
version (sin To = 1) is obtained on resonance when the 
pump power reaches the critical value. However, here, 
the critical value can be traversed without violating the 
validity of the equations. Full photon conversion is de- 
sirable in dynamical cooling: in that case, the higher 
frequency resonator will be emptied of photons, and the 
lower frequency resonator can be cooled to its ground 
state by pumping the intermediate frequency resonator 
(see lower panel of Fig. [2]). 



a ont [+ui] 




T aa t ac 




a in [+wi] 


_ c out [+us] _ 




tea ?*cc 




_ c in [+w 3 ] . 



where 



and 



(59) 
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2ip' 
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P >r 


2ip'* 




X^Xc 1 + 





^ = 1-1 



UJ 3 - UJ C 



-i(j> 



Vr Q r c 



(60) 

(61) 
(62) 



The reduced pump strength p' plays the same role here 
as p in the photon amplification case. Note that the 
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FIG. 4: Signal flow graph for a three-wave mixing device op- 
erating in conversion without photon gain, realized in case 
2 of Fig. [2] The incoming signal in port a (b) is reflected 
with amplitude r and transmitted with up-conversion (down- 
conversion) to port b (a) with amplitude (1 — r 2 ) 1 ^ 2 . 



C. Added Noise 

The number of output photons generated per mode in 
the amplification (case 1) is given by 



\r\ 2 K\ 



N 2 AC, 



(66) 



where Af m is the input photon spectral density given 
by Eq. (31) and we assume that there is no cross- 



correlations between the input fields a m and 6 11 



Assuming that the three-wave mixing device is in ther- 
mal equilibrium at temperature T <C hcoi^/kn and that 
the dominant noise entering the system at each port is 
zero-point fluctuations (j\f ln = 1/2), then in the 

limit of high gain \r\ 3> 1, the number of noise equivalent 
photons effectively feeding the system is 



M™=M°^/\r\ 2 ~l. (67) 

This means that the number of noise equivalent photons 
added by the device to the input is given by J\[ add = 
N™ — Af ln = 1/2. Hence, when operated as a non- 
degenerate amplifier with Gq 1, the device adds noise 
which is equivalent to at least half a photon at the sig- 
nal frequency to the input, in agreement with Caves 
theoreirPD. 

In contrast, in the conversion mode of operation, as- 
suming that there is no correlation between the input 
fields, the number of generated output photons per mode 
reads 



KT = \r\ 2 K n b + \t\ 2 ^ a - (68) 
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FIG. 5: Three-wave mixing element (see ellipse marked K 
in Fig. [TJ consisting of a loop of four nominally identical 
Josephson junctions threaded by a flux in the vicinity of half 
a flux quantum. Mutual inductances, not shown here, couple 
this circuit to inductances L a , Lt and L c of Fig. [l]via the 
inductances Lx, Ly and Lz respectively, which are much 
larger than the junction inductance Lj. The three currents 
Ix, Iy and Iz correspond to the three orthogonal modes of 
the structure. 



Therefore, in pure conversion where |r| = and \t\ = 1, 
when referring the noise back to the input, one gets noise 
equivalent photons 



■Nlq = A/" out / \t\ 2 = 1/2. 



(69) 



This means that, as a converter, the device is not required 
to add noise to the input since = W m . 



II. THREE- WAVE MIXING USING JRM 

The Josephson ring modulator is a device consisting 
of four Josephson junctions, each with critical current 
Io = 2iT7 forming a ring threaded by a flux $ = $0/2 
where <E>o is the flux quantum (see Fig. [5]) . The device has 
the symmetry of a Wheatstone bridge. There are thus 
three orthogonal electrical modes coupled to the junc- 
tions, corresponding to the currents Ix, Iy and Iz flow- 
ing in three external inductances Lx , Ly and Lz that are 
much larger than the junction inductance Lj = ip^E^ 1 , 
where ipo — h/2e is the reduced flux quantum. Each junc- 
tion j £ {a, P, 7, 6} is traversed by a current lj and at 
the working point (i.e. <E> = $o/2) its energy is, keeping 
terms up to order four in lj, given by 



Ej = -L e / f I* 
j 2 J 3 



1 r eff 

772 X 3 ' 
J 



24 



(70) 



where L e /^ = V%Lj and I' = Io/y/2. The currents ii 



the junctions are 


expressed by 
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(74) 



where J$ is the supercurrent induced in the ring by the 
externally applied flux $. The total energy of the ring 
is, keeping terms up to third order in the currents^, 



Ering — ^L j 



eff 



T 2 
L X 



Iy 



7 H 



1 L e / f U 
4 I'i 



We can express the currents as 

$a,&,c Ma,b,c 



l X,Y,Z 



L a ,b,c L x ,y,z L e JS ' 



Ixlylz- 
(75) 

(76) 



where M a ^. c are the mutual inductances between Lx,y,z 
and the oscillator inductances L a ^ iC . The non-linear co- 
efficient in the energy is, therefore, 



K = 



( £ 3") ! 



4<A) Ll ff L e h fS L e J f 



(77) 
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and we finally arrive at the result 



2 _ PaP b p c U}gUJ b Ld c 
93 ~ , eff 



Here the participation ratios are defined as 



-eff 



Pa,b,c 



eff 



and, at <I> = $ /2, 



U! 



eff 
J 



L 



128 Ej 



(78) 



(79) 



(80) 



The participation ratios are linked to the maximal num- 
ber of photons in each resonator, defined as those corre- 
sponding to an oscillation amplitude reaching a current 
of Iq in each junction of the ring modulator, 



Pa,b,c n a ,b.c 



T^a.b.i 



flUJ a . 



b.c 



(81) 



where the Ej' b ' c are of order Ej with factors accounting 
for the different participation of modes X, Y and Z in 
the current of each junction. Equations (78) and (81) 
are valid for all types of coupling between the Josephson 
ring modulator and signal/pump oscillators, which can 
be realized in practice by inductance sharing rather than 
by the mutual inductances discussed here. 

can also be rewritten in terms of the 
ating power in cavities a and b as 



Equation (81 



maximum circu 



pmax 
cav 



7a, b Ej 
Pa.b \f2 



(82) 



where we substituted Ej/y/2 as an upper bound for E a f b . 
The maximum number of photons in equation ( |81[ ) de- 
termine the maximum signal input power handled by the 
device 



a : b = £,la,bnUJa,bn afi . 



(83) 



We can now combine the notion of maximum power in 
resonator c compatible with weak non-linearity with that 
of a critical power for the onset of parametric oscillation 



given by Eq. (55): 



E 



- > nr - riF2 



p c fuv, 

arriving at the important relation 

PaPbQaQb > 3 



2 ' 



(84) 



(85) 



where 3 is a number of order unity depending on the 
exact implementation of the coupling between the ring 



Parameter 


Range 


UJ a ,b/2lV 


1 - 16 GHz 


Qa,b 


50 - 500 


Z a ,b 


10 - 150 n 


1c 


0.5 - 10 GHz 


Io 


0.5 - 10 fi A 


Ej 


10 - 230 K 


Pa,b,c 


0.01 - 0.5 


03/27T 


0.1 - 15 MHz 


fj max 


20 - 10 4 



TABLE II: Typical values for Josephson three-wave mixing 
devices. 



modulator and the oscillators. The quality factors of the 
resonators obey the well-known relation 



la,b 



u a ,b 

7a,& ' 



(86) 



Another maximum limit on the gain of the amplifier 
is set by the saturation of the device due to amplified 
zero-point fluctuations present at the input given by 



E, 



T ZPF 



a . b 



(87) 



Eqs. (81), (83) and (85) show that it is not possible to 
maximize simultaneously gain, bandwidth and dynamic 
range. 

In table II we enlist general bounds on the character- 
istic parameters of the three-wave mixing device, which 
are feasible with superconducting microwave circuits and 
standard Al-AlOx-Al junction fabrication technology. A 
few comments regarding the values listed in the table are 
in order. The frequency ranges of resonators a and b is 
mainly set by the center frequency of the system whose 
signal one needs to amplify or process. It is also im- 
portant that these frequencies are very small compared 
to the plasma frequency of the Josephson junction. The 
total quality factor range listed in the table (50 — 500) 
is suitable for practical devices. Quality factors in ex- 
cess of 500 can be easily achieved with superconducting 
resonators but, as seen from Eq. (58), higher the qual- 



ity factor, smaller the dynamical bandwidth of the de- 
vice. Quality factors lower than 50 on the other hand 
are not recommended either for a variety of reasons. For 
example, in the limit of very low Q the pump softens 
(becomes less stiff), and the dynamic range decreases as 
more quantum noise will be admitted by the device band- 
width and amplified "unintentionally" by the junctions. 
The characteristic impedance of the resonators Z a ^ is 
set by microwave engineering considerations as discussed 
in Sec. IV but, in general, this value varies around 50 
tt. The rate j c at which pump photons leave the circuit 
varies from one circuit design to the other as discussed 
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in Sec. IV and is limited by w c . This parameter also 
affects the maximum input power performance of the de- 
vice as explained in Sec. III. As to the values of Iq, on 
the one hand it is beneficial to work with large Josephson 
junctions in order to increase the processing capability 
of the device; on the other hand a critical current larger 
than 10 fiA adds complexity to the microwave design of 
the resonators and makes the fabrication process of the 
Josephson junction more involved. This might even re- 
quire switching to a different fabrication process such as 
Nb-AlOx-Nb trilayer junctions^ or nanobridgea^. The 
other parameters listed in the table, namely p a ,b,c, Ej, 
53> "'"fc X c their values depend, to a large extent, on the 
device parameters already discussed. 



III. LIMITATION OF DYNAMIC RANGE DUE 
TO PUMP DEPLETION 

In the last two sections, we were using results obtained 
by solving only the first two of the equations of motion 
Eqs. (26) under the restriction of the stiff pump ap- 



proximation. In this section, we extend our analysis and 
include the third equation describing the dynamics of the 
pump to calculate the pump depletion and its effect on 
the dynamic range of the device. For this purpose, we 
consider the average value of the third equation of motion 
for field c 

1 (c) = -zc c (c) - z 53 (ab) - ^ (c) + V% (c in (f)) . (88) 
In steady state and using RWA we obtain 

?<c(*))=V% <<?»(*)>■ (89) 



ig3(ab) + f(c(t))= y /^(S B (t)>. 

In the limit of vanishing input, the cross-correlation term 
(ab) is negligible and, therefore, 

2 



(c(t)) 



He 



(c in (i)) 



(90) 



The average number of photons in the c resonator in this 
case is, thus, 



lim fi c = -\(c in (t))\ 



(91) 



We now establish a self-consistent equation for n c , taking 
into account input signals of finite amplitude. We first 
evaluate the value of (a (t) 6(f)) in the frame rotating 
with the pump phase, 



(a(t)ft(t)) 



1 

2^ 



-foo /*+oo 



(a [u] b [a/]) e-^+^Mwdu/. (92) 



Using the field relations (see Appendix) 
^f a a M = a in [cj] + a out [uj] 



Jbb[uj] 



(93) 
(94) 



and the input-output relations given by Eq. (43), we 
obtain (transforming back into the time domain) 

leff (G) 



-ig 3 { a (t)b(t))=-^^-(c(t)) 
where, in the limit of large gains G> 1, 



(95) 



Jeff (G) 



2tt An" 



dc(A/TM+A/-rM)G(A w ) 

(96) 

denotes an effective decay rate of pump photons due to 
generation of entangled signal and idler photons. This 
last relation expresses, in another form, the Manley-Rowe 
relations^ that establish the equality between the num- 
ber of created signal photons by the amplifier to the num- 
ber of destroyed pump photons. It shows that even in the 
absence of any deterministic signal applied to the oscilla- 
tor a or b, pump photons are used to amplify zero-point 
fluctuations. Therefore, the pump tone always encoun- 
ters a dissipative load even when no signals are injected 
into the device. 

For a continuous wave (CW) input power sent at the 
center frequency of the a or b oscillator, or both, we have 



leff 



(Go,P in ) = ^G P 



477': 



(97) 



where P m = P" 1 + P™ is given in units of photon number 
per unit time and, in steady state, 



(G ,P ir 



47c 



( 7c + 7e// (G ,P ; ")r 



(98) 



As a finite input power is applied to the signal oscillators, 
oscillator c depopulates and, keeping the pump power 
constant, we get 



(G ,P h 



n c (G o ,P in = 0) 



1 



, G P'» 

G P in 
277' n 



(99) 



(100) 



On the other hand, from Eqs. (50) and (54), the left 
hand side is given by 



(G ,P h 



n c (G o ,P in = 0) 



VG-l 
\/G+l 



(101) 



where G denotes the gain in the presence of p n . In the 
large gain limit, if we fix the maximum decrease of gain 
due to pump depletion to be 



G_ 
Go 



> 1 



with e <C 1, then we obtain 



pir 

2nfl 



< eG, 



-3/2 



(102) 



(103) 



10 



which can also be rewritten as 



2n" 



GnP" 



> e 



(104) 



This relation shows that the ratio of the power of the 
pump tone to that of the signal at the output of the 
amplifier must always be much larger than the amplitude 
gain, in order for the linearity of the amplifier not to be 
compromised by pump depletion effects. 

In Fig. [6] we plot a calculated response of the signal 
output power P out versus the signal input power P; n for 
a typical three- wave mixing device. The device parame- 
ters employed in the calculation and listed in the figure 
caption are practical values yielding a maximum input 
power, which is limited by the effect of pump depletion. 
The different blue curves are obtained by solving Eq. ( 98 ) 
for G and using the input-output relation P out = Gp n , 
where P n expressed in units of power is taken as the 
independent variable and Go is treated as a parameter. 
Note that in solving Eq. Q, equations ([97]), ([91]) , ( [50] ) 
and (54 1 are used. When drawn on logarithmic scale, 



the device gain translates into a vertical offset (arrow in- 
dicating Go) off the P ou t = Pj n line, indicated in red. 
The dashed black vertical line corresponds to a signal 
input power of 1 photon at the signal frequency per in- 
verse dynamical bandwidth of the device at Go = 20 
dB. The dashed green line corresponds to the maximum 
gain set by the amplified zero-point fluctuations given by 
Eq. (87), while the cyan line corresponds to the maxi- 
mum circulating power in the cavity given by Eq. ( |82| . 
Furthermore, the maximum bound p max indicated by 
the solid magenta line corresponds to P£^ x = GoP i ' n lax , 



where P™ ax = P 6 ma> 



/Go 



3/2 



and P h max = P™ ax . As can 



be seen in the figure the predicted power scaling due to 
pump depletion effect, expressed in relation (103), follows 



the calculated response quite well. Finally, it is straight- 
forward to see that the usable region in the parameter 
space of the device with respect to gain, bandwidth and 
maximum input power lies within the boundaries of the 
fictitious triangle ABC indicated in the figure which is 
formed by the intersection of the magenta, black and red 
lines. 



IV. THE JOSEPHSON PARAMETRIC 
CONVERTER 
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FIG. 6: (Color online). A calculated response of the signal 
output power P out versus the signal input power P; n of a typ- 
ical three-wave mixing device which exhibits a pump deple- 
tion effect. The different blue curves correspond to different 
Go setpoints. The definition of the other lines in the figure 
is given in the text. The parameters used in the calcula- 
tion are: Lu a /2iv = 7 GHz, u b /2ir = 8 GHz, uj c /2tt = 15 GHz, 
-y a /2n = j b /2iv = 50 MHz, 7 c /2tt = 0.6 GHz, Q a = 140, 
Q b = 160, p a = Pb = 0.03, Pc = 0.02, Jo = 1/*A, 



E" 



Ej/s/2 



16.3 K, P 



lph 



-128 dBm, G? 



35 



dB, P c ^ x = P 6 max = -86 dBm and g 3 /2n = 0.7 MHz. 



A. Microstrip Resonator JPC (MRJ) 



Mode b a 
O 



50 Q 




o o 

Mode c Mode a 



We discuss here three different realizations of the 
Josephson parametric converter (JPC), which constitutes 
a fully non-degenerate three-wave mixing device capable 
of amplification and conversion as discussed in the pre- 
vious sections. The three schemes differ in the resonator 
circuit design and the coupling between the feedline and 
the resonator. 



FIG. 7: Circuit model of the Microstrip Resonator JPC 
(MRJ). 



The MRJ comprises two superconducting microstrip 
resonators which intersect at a JRM at the center as 
shown in the circuit model of the device in Fig. [7j The 



11 



resonance frequencies of the MRJ are determined by the 
lengths of the microstrips l a — A a /2 and 4 ~ A&/2 and 
the Josephson inductance of the JRM, where A a , A& are 
the wavelengths of the fundamental resonances at oj a and 
uji,- It is worth mentioning that in addition to the differ- 
ential modes a and 6, this configuration of two coupled 
resonators also supports a common (even) mode. The 
angular frequency oj e at which this even mode resonates 
lies between (uj b + cu a ) /2 and u>b (where Wb > w a ). The 
characteristic impedance of the resonators in the MRJ 
model is designed to be 50 £1 to ensure optimal coupling 
to the feedlines. Figure [8] exhibits an optical image of 
a typical MRJ device. The resonators are usually made 
of Al or Nb over sapphire or high-resistivity silicon and 
are coupled to the (transmission-line) feedlines using gap 
capacitors. The main role of these coupling capacitors is 
to set the external quality factor of the resonators. For a 
large bandwidth device operating in the 6 — 10 GHz band, 
the external Q of the resonators is typically in the range 
60 — 100. In all JPC designs discussed here the total Q 
essentially coincides with the external Q, since the inter- 
nal losses of the resonators are less than 10~ 4 . Signals at 
oji and u)2 , which lie within the bandwidths of resonators 
a and 6, are fed into the JPC through the delta port of 
a 180 degree hybrid, whereas the pump drive applied at 
lu 3 = uji +u; 2 , f° r amplification, is a non- resonant tone 
and is injected into the device through the sigma port of 
the hybrid (Fig. 0. The main advantage of the MRJ 
is that it is easy to design and fabricate. On the other 
hand, the main disadvantages are: (1) the area of the de- 
vice can be relatively large depending on the frequencies 
of interest, (2) the characteristic impedance of the de- 
vice is limited to around 50 (3) the pump can be less 
stiff than the designs discussed below. The latter is due 
to the fact that the transmission-line resonators support 
higher resonance modes such as 2uj a and 2u b with finite 
Q, which can be relatively close to the pump angular 
frequency lo 3 . 



B. Compact Resonator JPC (CRJ) 

In order to mitigate some of the drawbacks of the MRJ, 
we developed a new JPC design based on compact res- 
onators known as CRJ. The circuit model of the CRJ, 
shown in Fig. [9] consists of four equal capacitors denoted 
as C and two pairs of linear inductors connected in se- 
ries with the JRM whose total inductance is L a and L b 
respectively. Using symmetry considerations one can ver- 
ify that this circuit has three eigenmodes. Two differen- 
tial eigenmodes which resonate at bare angular frequen- 

cies u a = l/^L a + L e / f ^C, io b = 1/J + L°/ f> ) C, 

where Ly is the equivalent Josephson inductance of the 
JRM biased at half a flux quantum, and an even eigen- 
mode which resonates at a lower bare ang ular frequency 
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FIG. 8: (Color online). Optical microscope image of a mi- 
crostrip resonator JPC (MRJ). The resonators denoted a and 
b are half-wave microstrip resonators which intersect at a 
JRM. A zoomed-in view of the JRM, which consists of four 
Josephson junctions arranged in Wheatstone bridge configu- 
ration, is shown on the right. The MRJ is coupled to 50 Q 
feedlines via gap capacitors. 
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C. Figure 
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shows an opti- 



FIG. 9: Circuit model of the Compact Resonator JPC (CRJ). 



cal image of a typical compact JPC. The resonators of the 
device are made of Nb deposited over sapphire substrate. 
They are fabricated using a standard photolithography 
step and RIE etching. The JRM at the center of the de- 
vice is made of Aluminum. It is fabricated using e-beam 
lithography, and angle shadow evaporation. As can be 
seen in the figure, the capacitance elements (including 
the coupling capacitors) of the device are implemented 
using interdigitated capacitors, whereas the inductive el- 
ements are realized using long narrow superconducting 
lines. Unlike the microstrip resonator JPC, the compact 
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resonator JPC does not have higher harmonic resonances. 
The next closest resonance of this structure resides above 
4w a , therefore the pump applied at uj a + ojf, can be con- 
sidered stiff to a very good approximation. Other advan- 
tages of this realization are: (1) small size, with dimen- 
sions much smaller than the wavelengths corresponding 
to the resonance frequencies, (2) no requirement of a defi- 
nite ground plane, unlike the MRJ, (3) greater flexibility 
in engineering the characteristic impedance of the res- 
onators higher or lower than 50 f2, (4) higher internal 
quality factor resonators than the microstrip design. On 
the other hand, the main disadvantages of this design 
are: (1) the narrow lines and the interdigitated capaci- 
tors (as well as the lines connecting them) have parasitic 
capacitances and parasitic inductances associated with 
them, therefore scaling these devices to match a certain 
frequency or certain characteristic impedance requires us- 
ing a microwave simulation tool, (2) there is a limit to 
how big the capacitance can be using the interdigitated 
configuration (values above 0.5 pF is difficult to achieve), 
therefore engineering characteristic impedances below 30 
is not quite feasible with this design. 




FIG. 10: (Color online). Optical microscope image of a com- 
pact resonator JPC (CRJ). The device consists of four equal 
interdigitated capacitors denoted C and two inductive ele- 
ments denoted L a and Lb which are realized using narrow 
superconducting lines of different lengths. The JRM of the 
device resides at the intersection of the two lines. An optical 
image of the JRM is shown in the inset. The CRJ is cou- 
pled to 50 microstrip feedlines via interdigitated capacitors 
denoted C c . 



C. Shunted JPC (SJ) 

In this subsection we discuss a third promising design 
called the capacitively and inductively shunted JPC (SJ) 
which is still a work in progress in our lab. In this version 




A Mode a 



2 Mode c 



FIG. 11: Circuit model of the Shunted JPC (SJ). 



of the JPC, the capacitive elements are parallel plate ca- 
pacitors and the inductive elements are mainly Joseph- 
son junctions. A schematic circuit model of the SJ is 
drawn in Fig. 



11 It is straightforward to show that the 



SJ model has two differential eigenmodes with angular 

resonance frequencies ui a = 1/ J L'jC a , cot, = 1/ J L'jCb, 

where L j corresponds to the equivalent inductance of 
the JRM shunted by linear inductor^!, as shown in Fig. 



11 The main purpose of these shunting inductors is to 
eliminate the hysteretic flux response of the JRM and 
extend the frequency tunability of the device beyond the 
bandwidth limit of the resonators. Such frequency tun- 
ability is achieved by varying the flux threading the loop 
which, in turn, varies Lj. Note that the addition of these 
shunting inductors can be employed in other realizations 
of the JPC also, such as the MRJ, as shown in RefpDand 
the CRJ. It is important to emphasize, however, that the 
main difference between the SJ and the CRJ or MRJ 
schemes is that the shunted JRM in the SJ design is the 
only inductive element in the circuit that forms an in- 
tegral part of the resonators a and b. Thus, the larger 
lumped capacitors employed in the SJ design play a cru- 
cial role in keeping the resonance frequencies of the device 
below 10 GHz. 

Similar to the Josephson bifurcation amplifier (JBA) 
implementation^, the plate capacitors in the SJ design 
can be made of Nb electrodes separated by a thin SiN 
dielectric layer. Using plate capacitors in this realization 
has two advantages: (1) the plate capacitors can be made 
very large, i.e. their capacitance can vary in the range 
1 — 40 pF, (2) they are easy to design as their capacitance 
scales linearly with the electrode area. Furthermore, due 
to the lumped nature of the capacitive and inductive el- 
ements in the SJ design and the fact that the capacitors 
can be large, the SJ has three important advantages over 
the previous designs: (1) the characteristic impedance of 
the resonators can be of the order of a few ohms, which 
yields an improved coupling between the resonators and 
the JRM, (2) due to the impedance mismatch between 
the characteristic impedance of the resonators and the 
50 £1 feedlines, the coupling capacitors are unnecessary 
to achieve low external Q and the feedlines can be con- 
nected directly to the resonators, (3) the maximum input 
power of the amplifier can be increased by increasing the 
critical current of the JRM junctions while keeping the 
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resonance frequencies fixed by enlarging the capacitors. 



V. EXPERIMENTAL RESULTS 



the working point Ej = Ioipo/^/2, where Iq is evaluated 
using dc resistance measurement of the junctions. 

Using Eqs. ( 50 ) and ( 54 ) for the case of maximum gain 
Gj? ax yields 



The set of JPC parameters which can be directly mea- 
sured in an experiment are: the angular resonance fre- 
quencies of the resonators a and b u) a , u>b, the inverse of 
residence times of photons at resonance j a , 7b, the par- 
ticipation ratios p a , pb, the maximum input power which 
the device can handle with no applied pump tone P™ ax , 
P b max , and the maximum measured gain at vanishing in- 
put power G nax . 

One way to find p a , pb is by measuring uj a , cj& as a func- 
tion of applied magnetic flux threading the JRM loop. 
To establish this relation, we model the resonators near 
resonance as an LC oscillator with effective inductance 
Laj, and effective capacitance C a .b- In this model, the 
bare angular resonance frequencies of the device (with 
the junctions) uj a , ujb, can be written as 



Ua,b (<P) = 



y/C a ,b (L a ,b + Lj (ip))' 



(105) 



where Lj (ip) is the effective Josephson inductance of the 
JRM given by 



Lj(<p) 



L, 



cos 



(!) 



(106) 



with tp — 27r$/$o. By calculating the derivative of 
u a ,b (v) with respect to the reduced flux ip, one gets 



1 du a ,b 

Wa,i dip 



L.J (V) 



(107) 



1 ftp\ 
- tan — . 

8 \aJ (L a>b + Lj (ip)) 
|tan(|)p , 6 fo>). (108) 



Hence, at the device working point $ = $o/2 (<ys = 7r), 
p a ,b reads 



Pa,b = 



1 duj a: i 
u a ,b dip 



(109) 



1^ — 7T 



Furthermore, using Eq. (81) and the measured values 
u., x Dmax^ one can m f er ^h e Josephson energy E a j 



which is available for amplification 



Ej = Pa,b 



pmax 
r a,b 

7a,b 



(110) 



It is important to mention that, in our experiments, 
we find that this value is lower by about one order of 
magnitude than the Josephson energy of the junctions at 



2- _ 7a lb \[G 

93 n c,p-yl — 



max 




(111) 



which in the limit of high gains gives an upper bound on 
the product <?|n CiP _s.i 



< 



7 alb 

4 



(112) 



Here n c p _>i is the number of pump photons in the device 



at G^ ax . In Figs. (12), (13 1, (14) we plot on logarithmic 
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FIG. 12: (Color online). Output power P out measurement 
of a CRJ amplifier (device A) as a function of input power 
Pin measured at u>b- The data curves plotted in blue corre- 
spond to different Go setpoints obtained for different pump 
powers. The red line corresponds to dB (unity gain) where 
Pout = Pin- The dashed black vertical line indicates the in- 
put power of 1 photon at the signal frequency per inverse 
dynamical bandwidth of the device at Go = 20 dB. The top 
horizontal line labelled P^v X corresponds to the maximum 
circulating power in the resonator cavity given by Eq. ( 82 1 . 



The green line corresponds to an upper limit on the device 
gain set by the saturation of the amplifier due to zero-point 
fluctuations given by Eq. ( |87[ ). The dashed magenta line is 
a theoretical prediction for Po^t X ) which corresponds to the 
maximum circulating power in the device given by Eq. ( 114 1 



The measured and calculated parameters of this device (A) 
are listed in table III. 



scale the output power P out of three different JPCs with 
different characteristics as a function of input power P ln . 
For simplicity, we refer to the three devices as A, B and 
C respectively. The parameters of the three devices are 
listed in table III. The data curves plotted in blue are 
measured at resonance and satisfy the relation 



-Pout — G (P; n , Go) Pi, 



(113) 
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where G (Pin, Go) is the amplifier gain. This depends 
on Pi n and Go, the device gain for P; n = which is set 
by the applied pump power. In this measurement, we 
apply a fixed pump power and vary P n treating Go as a 
parameter. In log units, the device gain translates into 
a vertical offset from the dB baseline (red line) which 
corresponds to P out = p n - 

As expected, the devices maintain an almost constant 
gain Go as a function of P n before they saturate and 
their gain drops for elevated input powers. However, as 
can be seen in Figs. (12), (13), (14), the three devices 



exhibit qualitatively different behaviors in the vicinity of 
their maximum input power, which correspond to differ- 
ent saturation mechanisms taking place in the device as 
will be discussed shortly. Note that the order in which 
the different results are presented in this section does 
not depend on the specific implementation of the device 
(see Sec. IV) but rather on the saturation mechanism 
involved in each case. 

In Fig. 12 device A exhibits almost a plateau in P out 
as it reaches its maximum input power for different Go 
setpoints. This result can be explained by assuming a 
stiff pump for which Eq. (83) applies. By employing 



P h max , measured with no applied pump tone, we plot the 
dashed magenta line labelled P max which corresponds to 

cr = pr x - (H4) 

The dashed black vertical line indicates the input power 
of 1 photon at the signal frequency per inverse dynamical 
bandwidth of the device at Go = 20 dB. In practice, as 
we discuss in Sec. VI, the usable region in the parameter 
space of the device with respect to gain, bandwidth and 
maximum input power lies within the boundaries of the 
fictitious triangle formed by the magenta, red and black 
lines. Furthermore, in Figs. (12), (13), (14) we plot two 



fundamental limits on the maximum gain G^pp (green 
line) which corresponds to saturation of the device due 
to amplified zero-point fluctuations and the maximum 
circulating power P"™ K (cyan line), given by Eq. (87) 
and Eq. ( 82 ) respectively. 



The fact that these lines lie considerably above the 
experimental data in Figs. (12), (13), (14), suggests 



that the energy threshold, at which nonlinear effects in 
these devices become significant, is much lower than the 
Josephson energy of the junctions, i.e. Ey b <C Ej. 

In contrast to Fig. 12 the data curves shown in Fig. 
[T3| for device B, exhibit a gradual decrease in the gain 
in the vicinity of the maximum input power which can 
be explained in terms of pump depletion effect discussed 
in Sec. III. The maximum bound P max indicated by 
the solid magenta line corresponds to P™^ x = G P^ ax , 
where in this case P i '" ax 
is given by 



satisfies the inequality 103 and 



pu 



G 



3/2 



(115) 




FIG. 13: (Color online). Output power P ou t measurement of 
a MRJ amplifier (device B) as a function of input power P n 
measured at uj a - The data curves plotted in blue correspond 
to different Go setpoints obtained for different pump powers. 
The red line corresponds to dB (unity gain) where P out = 
Pin- The dashed black vertical line indicates the input power 
of 1 photon at the signal frequency per inverse dynamical 
bandwidth of the device at Go = 20 dB. The top horizontal 
line labelled P c ™v X corresponds to the maximum circulating 
power in the resonator cavity given by Eq. ( |82[ ). The green 
line corresponds to an upper limit on the device gain set by the 
saturation of the amplifier due to zero-point fluctuations given 
by Eq. (87 1. The solid magenta line is a theoretical prediction 
for P„ ax of the device and the cor resp onding P^ x due to 
pump depletion effect given by Eq. ( 115 1. The measured and 



calculated parameters of this device (B) are listed in table III 



device C exhibit an abrupt drop in the device gain in the 
vicinity of P™ ax of the device, which indicates that the 
device enters an unstable regime at elevated input pow- 
ers. As can be seen in this case the solid magenta line 



which satisfies P^ ax 



G P™ ax , where P™ ax is given 
by Eq. jll5|) — lies above the experimental data. This 
suggests that the maximum input power in this sample, 
which displays a steeper power scaling than Eq. (115), 



is mainly limited by nonlinear effects arising from higher 
order terms in the Hamiltonian of the system and can- 
not be attributed to a pump depletion effect alone. It 
is worthwhile noting that a similar power scaling for the 
maximum input power has been observed as well for an 
MRJ amplifier in Ref. H2 

To understand which properties are responsible for the 
different gain behaviors exhibited by devices A, B and C, 
we point out a few important distinctions in their design 
(respective parameters are listed in table III) . The data in 
Fig. 12 (device A) and Fig. 14 (device C) is measured on 



On the other hand the data curves shown in Fig. 14 for 



JPC devices realized using the CRJ configuration which 
yields, in general, a stiff pump response as explained in 
Sec. IV (B). However, the main two differences between 
devices A and C are: (1) device A has a narrower band- 
width as compared to C (70 MHz vs. 142 MHz) and (2) 
the JRM junctions in A have a smaller Iq compared to 
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FIG. 14: (Color online). Output power P ou t measurement of 
another CRJ amplifier (device C, with different parameters 
from device A) as a function of input power P n measured at 
u)b- The data curves plotted in blue corresponds to different 
Go setpoints obtained for different pump powers. The data 
curves of this device exhibit abrupt drop in the gain in the 
vicinity of the maximum input powers which suggests that 
the device enters an unstable regime at elevated input powers. 
The red line corresponds to dB (unity gain) where P ou t = 
Pi n . The dashed black vertical line indicates the input power 
of 1 photon at the signal frequency per inverse dynamical 
bandwidth of the device at Go = 20 dB. The top horizontal 
line labelled P^v" corresponds to the maximum circulating 
power in the resonator cavity given by Eq. (821. The green 



line correspond to an upper limit on the device gain set by the 
saturation of the amplifier due to zero-point fluctuations given 
by Eq. (87 1. The solid magenta line is a theoretical prediction 
for Pi™ a of the device and the corresponding P^ ax due to 
pump depletion effect given by Eq. 1 115 1. The measured and 



calculated parameters of this device (C) are listed in table III 



those in C (2 /iA vs. 4 /iA). The relatively large band- 
width of device C leads to a larger dynamical bandwidth 
14 MHz at G = 16 dB, as opposed to 10 MHz achieved 
in device A for the same gain, and also yields (with the 
larger 1$ of device C) higher P^ x values. However, the 
large bandwidth translates into a lower pQ product for 
C as compared to A, thus making it more susceptible to 
parametric oscillation (at high gains or high input pow- 
ers) as implied by inequality ( 85 1 . 



Device B, on the other hand, exhibits a pump depletion 
effect as shown in Fig 



13 



This can be attributed to its 
MRJ configuration, which, in general, exhibits a less stiff 
pump response than the CRJ, due to the presence of high 
order modes as explained in Sec. IV (A). Furthermore, 
as opposed to the MRJ amplifier in Ref. [12] with an 
idler frequency of 6.4 GHz, device B has a higher idler 
frequency of 15 GHz which leads to a higher pQ product. 



Parameter \ Device 


A 


B 


C 


Design 


CRJ 


MRJ 


CRJ 


u a /2iT (GHz) 


6.576 


8.436 


7.051 


OJ6/27T (GHz) 


6.873 


15.087 


7.673 


cu 3 /2ty (GHz) 


13.449 


23.523 


14.724 


7a /27T (MHz) 


69 


116 


79 


7t,/27r (MHz) 


71 


250 ± 25 


142 


Qa i Qb 


94,96 


73,60 


89,54 


Pa,Pb 


0.02 


0.03,0.05 


0.03 


PaPbQaQb 


8.1 


6.6 


4.3 


i (M) 


2 


3 


4 


P- V ax (dBm) 


-82 


-77 


-76 


P a m b ax (dBm) 


-97 


-89 


-87 


Pi ph (dBm) 


-127 


-125 


-123 


GS (dB) 


38 


39 


40 


G max (dB) 


22 


20 


16 


gsnJ^^TT (MHz) 


33 


77 


45 



TABLE III: Parameters of JPCs A, B and C. Precision is last 
significant digit unless indicated otherwise. 



VI. REQUIREMENTS FOR QUBIT READOUT 

One of the leading architectures which is used to ma- 
nipulate and readout the state of s uperc onducting qubits 
such as transmons and fluxoniuma^ESl is circuit Quan- 
tum Electrodynamics (cQED). In such a system a quan- 
tum non-demolition measurement of the qubit state can 
be performed using dispersive readout in which the fre- 
quencies of the qubit and the cavity are detuned. As a 
result, the qubit and the cavity interact via exchanging 
virtual microwave photons^ and the qubit state gets en- 
coded in the output microwave field of the cavity. How- 
ever, since the energy of microwave photons is very small, 
the detection of single photons is difficult especially con- 
sidering the fact that state-of-the-art cryogenic ampli- 
fiers (i.e. high electron mobility transistor (HEMT) 35 ) 
following the cQED setup add noise to the input signal, 
equivalent to about 20 — 40 photons at the signal fre- 
quency. Therefore, adding a quantum-limited amplifier 
in series between the cQED sample and the HEMT am- 
plifier can substantially decrease the noise temperature 
of t he sy stem and enable real-time tracking of the qubit 
stat d 9 * 10 l The desired requirements of a Josephson para- 
metric amplifier for such high-fidelity qubit readout can 
be summarized as follows: 

• A center frequency in the range 5—12 GHz which is 
widely used in readout cavities of superconducting 
qubits. 

• A large power gain on the order of 20 dB in order 
to beat the noise of the following amplifier, i.e. the 
HEMT. 

• A minimum added noise to the signal, equivalent to 
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a half input photon at the signal frequency Tn = 
hu) a .b/2kB when operated in the phase preserving 
modeP. 

• A large dynamical bandwidth of the order of 10 
MHz, which corresponds to a signal processing time 
of less than 100 ns and matches the bandwidth of 
most readout cavities. 

• A maximum input power of a few photons per in- 
verse dynamical bandwidth of the device at the 
highest gain. Such requirement is essential in quan- 
tum non-demolition readout schemes which employ 
of the order of a photon on average^. 

• A tunable bandwidth of more than 100 MHz so that 
the center frequency of the amplifier can match the 
frequency of the readout signal. Recently, Roch et 

have achieved a tunable bandwidth of more 
than 500 MHz in a MRJ device by shunting the 
Josephson junctions of the JRM with linear induc- 
tors realized using superconducting wires [see Fig. 
Similar results were obtained by our group in a 



11 



MRJ device by utilizing large Josephson junctions 
instead of superconducting whea^S. 

• Minimal out-of-band back-action to avoid qubit re- 
laxation. 

In table IV we enlist the parameters achievable with a 
JPC which show its viability as a low-noise preamplifier 
for qubit measurements. 



Property 


Requirement 


Achieved to date^ 


0Ja,b/2n 


5 - 12 GHz 


6.4 & 8.1 GHz 


T N ® 8 GHz 


0.2 K 


0.4 K 


G 


20 <1B 


21 dB 


B 


10 MHz 


11 MHz 


Tunable BW 


100 MHz 


60 MHz 


p 

J max 


1 photon 


3 photons @ 20 dB 


OB back-action 


Negligible 


None measurable 



TABLE IV: Preamplifier requirements and JPC merits 
achieved to date (OB=out-of-band). 

The last question which we would like to address in 
this paper is whether there exists a set of technologically 
feasible parameters for which the JPC can be optimized 
with respect to dynamic range while maintaining a gain 
of 20 dB and a reasonable dynamical bandwidth of more 
than 2 MHz. In order to provide a quantitative answer 
we choose a signal frequency of 12 GHz, which is a good 
choice for readout frequency for qubits as it is higher 
than most qubit frequencies. We set an ambitious goal 
for the processing capability of the JPC of about 100 
input photons at the signal frequency with gain 20 dB. 
To that end, we choose to perform the optimization pro- 
cess for the CRJ configuration which inherently yields 
large 7 C values and also allows variation of the external 
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FIG. 15: (Color online). An optimized JPC response drawn 
in blue for large maximum input power in excess of 100 input 
photons at 12 GHz per inverse dynamical bandwidth of the 
device at 20 dB. The definition of the other lines shown in 
the figure is similar to Fig. ( |12[ |. The parameters used in the 
calculation are: u) a /2n = 11 GHz, u)b/2n = 12 GHz, uj c /2tt = 
23 GHz, Z a = 36 fi, Z b = 33 Q, L a = 0.51 nH, L b = 0.42 
nH, C a = C b = 0.4 pF, C Ca = 31 fF, C c „ = 28 fF, Ja /2ir = 
7b /27r = 44.2 MHz, 7c /2tt = 3 GHz, Q a = 249, Q b = 271, 
Qc = 8, p a = 0.028, p b = 0.034, p c = 0.02, I = 30 n A, 
Ej/V2 = 490 K, E a j' b = 49 K, P lph = -126 dBm, G^ = 48 



dB 



ria.b 

•jj = <±vr\, .riph 
= —86.3 dBm and gs/2n 



■■ 0.6 MHz. 



quality factor of the resonators more easily than the SJ 
scheme. We also choose a resonance frequency for mode a 
of LOa/2-K = 11 GHz, a relatively high critical current Iq = 
30 /iA and limit ourselves to capacitance values below or 
equal to 0.4 pF. The advantage of working with large Iq 
for the purpose of large dynamic range is that it increases 
Ej and lowers g%. However, such large Jo yields very low 
Ifjf = t/2ipo/Iq = 15 pH which requires coupling to rel- 
atively low impedance resonators while maintaining par- 
ticipation ratios of a few percent. The next challenge in 
the optimization is to increase the pQ product of the de- 
vice which promotes high gains by increasing the quality 
factor of the resonators. Nevertheless, care must be taken 
not to increase the quality factors beyond what is strictly 
necessary for two reasons (i) high Q resonators limit the 
dynamical bandwidth of the device as can be seen in Eq. 
(58), (ii) high Q resonators increase the pump depletion 
effect and, in turn, lower the dynamic range of the de- 
vice. In Fig. [l5]we plot the calculated response of such an 
optimized JPC which takes into account the above con- 
siderations and limitations. As can be seen in the figure 
the optimized device exhibits, for the chosen set of pa- 
rameters, a maximum input power of about 100 photons 
at the signal frequency per inverse dynamical bandwidth 
of the device B/2n = 4.4 MHz at 20 dB of gain. The 
device parameters which are used in the calculation are 
listed in the figure caption. It is important to note that 
in the calculation of the expected response, which is in- 
dicated by the blue curves for different values of Go, we 
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assumed an available Josephson energy 10 times smaller 
than Ej/^/2 of the junctions, in agreement with experi- 
mental conditions. Finally, we verify that the set of pa- 
rameters of the optimized device satisfy the inequalities 
iir = 3.7 • 10 3 > n?° = 1.3 • 10 3 > n 2 ° dB = 10 3 . 



VII. CONCLUSION 

We have addressed in this paper a new type of quan- 
tum signal processing device based on Josephson tunnel 
junctions. In contrast with the devices based on SQUIDS 
and driven non-linear Josephson oscillators, it performs 
a fully non-degenerate three-wave mixing in which the 
modes of the signal, pump and idler are separate both 
spatially and temporally. The heart of the device con- 
sists of a ring modulator constructed from four Josephson 
junctions arranged in a loop. Both quantum-limited am- 
plification and noiseless frequency conversion are possible 
with this device, and the characteristics of these analog 
signal processing operations are entirely calculable ana- 
lytically. We have established the limitations preventing 
the simultaneous maximization of photon number gain, 
bandwidth and dynamic range. Nevertheless, we have 
shown that a device satisfying all the requirements of 
superconducting qubit readout is realizable with present 
day technology. 
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FIG. 16: The damping of a circuit by a resistance R can take 
place in a parallel or series way, depending on whether the 
resistance is placed across a branch or in series with it. The 
Nyquist model represents the resistance by a transmission line 
with characteristic impedance Z c = R. The noise source as- 
sociated with the resistance (fluctuation-dissipation theorem) 
is a parallel current source in the parallel case and a series 
voltage source in the series case. The noise source is replaced 
in the Nyquist model by incoming thermal radiation whose 
amplitude A' u is the square root of the power flux of the radi- 
ation (A ln should not be associated to a vector potential and 
is rather like the square root of the length of the Poynting 
vector) . 
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Appendix: Quantum signals propagating along a 
transmission line and input-output formalism 

This appendix treats quantum-mechanically the damp- 
ing of a circuit by a resistance modelled as a semi-infinite 
transmission line, as shown in Fig. |16[ It borrows heavily 
from the book by Gardiner and Zollepni but uses slightly 
different notations that are adapted to the specificities 
of our Josephson circuits. We first describe an infinite 
transmission line extending from x — — oo to x = +oo. 
Later, we will cut the line at x — and replace the left 
portion by two terminals of the circuit. 



Infinite transmission line 

The capacitance and inductance per unit length of the 
line are Ct and Lt, respectively. The equations obeyed 



by the current / along and the voltage V across the line 
are 



(116) 
(117) 



in which, for the moment, we treat the fields classically. 
The characteristic impedance and propagation velocity 
are given by 



Z r = 



Li 

cV 
i 

LtCe 



(118) 
(119) 



In order to solve Eqs. (116) and (117), we introduce 
two new fields: the left-moving and right-moving wave 
amplitudes, 



(x,t) 
A^ (x, t) 



1 

i 

M 



V{x,t) + ^Z c I(x,t) 
V(x,t)-y/ZJ (x,t) 



(120) 
(121) 



which have the advantage of treating currents and volt- 
age on the same footing (note that these amplitudes are 
not directly related to the vector potential) . The dimen- 
sion of these fields is [watt] 1 / 2 and they are normalized 
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such that the total power P traversing, in the forward 
direction, a section of the line at position x and time t is 
given by 



P(x,t) = [A-+ {x,t)f-[A^ (x,t)} 2 . 



(122) 



The quantity P here plays the role of the Poynting vector 
in full 3D electrodynamics. Each of the terms at the 
right hand side of the last equation is thus the separate 
contribution of the corresponding wave to the total power 
flow. 

When solving Eqs. ( fTl6pTf| ), we find 



d lT ± , \ 1 d tT± , 
ox v p at 



(123) 



This relation means that A*^ does not depend separately 
on x or t but a combination of both and thus: 



= A^ (x — v p t, t = 0) . 



(x,t) =A^[x = 0,t 



A*- (x,t) = A*- [x = 0,t + 



A*~ (x + v p t, t = 0). 

(124) 



The properties of the wave amplitude can be summarized 
by writing 



A^ (x, t) = Aq (t) 



t = t H X, 



e*~ =Tl- 



(125) 
(126) 
(127) 



Note that the detailed definition of the retardation t de- 
pends on the wave direction. We now turn to the energy 
density U (x,t), related to P by the local energy conser- 
vation law 



dU _ _dP 
dt dx ' 



(128) 



Combining Eqs. (122) and (123), we get 
dU (x,t) 



dt 
2 



• v 

= ^{[A- {x,t)f + [A^ {x,t)f} . (129) 
The total energy of the line at time t is, thus^, 

\[A^ (x,t)] 2 + [A<- (x,t)] 2 }dx. (130) 



d d 
A (x, t) —A^ (x, t) + A*- (x, t) —A^ (x, t) 



H = 



+oo 



When H in Eq. (130) is considered as a functional of 



dynamical field variables A~^ and A*~ , the equation of 



motion Eq. ( 123 1 can be recovered from Hamilton's equa- 

(131) 



tion of motion as 
8 



dt 



A- (x,t) = - {H, A*~ (x,t)} 



P.B. 



on imposing the Poisson bracket 

-> 19 
{A- {x 1 ,t 1 ),A^- (x 2 ,t 2 )} pB = - _ 6(tx-t 2 ). 

(132) 

Therefore, from the classical-quantum correspondence in- 
volving the replacement of Poisson brackets by commu- 
tators, we find that the quantum operator version A*~ of 
the fields satisfy the commutation relation 



A^ (xi,ti) ,A^ (x 2 ,t 2 ) 







2 d{ Tl -T 2 ) 



S (n - t 2 ) , 
(133) 

which is analogous to the commutation relation between 
the electric and magnetic field in 3-D quantum electrody- 
namics. Note that the fields are Hermitian at this stage. 
Introducing the Fourier transform, 



A* M = 



i 



'2lt 



+ GO 



A^ (x = 0, r) e iwT dr, (134) 



where the Fourier components (which are now non- 
hermitian operators) satisfy 



A** H f = A* [-«] , 



we can also write the Hamiltonian as 

«+oo 



(135) 



E 



A a [u] A a [-u] duo. 



(136) 



The field operators in the frequency domain satisfy 



A^ [wi] , A^ [uj 2 ] = - (u)i — oj 2 ) S (wi + lo 2 ) . (137) 



We now introduce the usual quantum field annihilation 
operators 



a M = 



i- [u] 



= a — lo 



t 



A 

They satisfy the commutation relations 



u \ <- r it 

a' \LU\ = — = a 



(138) 
(139) 



[a^ [wi] [u 2 ]\ = sgn 



It is useful to note that since 



uii - uj 2 



5 [ui\ + uj 2 ) . 



a \uj = a —a; 



(140) 
(141) 
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Eq. ( 140 1 exhaustively describes all possible commutator 
cases. 

In the thermal state of the line, at arbitrary tempera- 
ture (including T = 0), 



(a^ [wi] [uj 2 ]) = S a *± a *± 
where 



2 



5 (wi + uj 2 ) , 

(142) 

S ^M=sgn(w)JVT(w). (143) 

When cj is strictly positive Nt (w) is the number of avail- 
able photons per unit bandwidth per unit time travelling 
on the line in a given direction around frequency w 



N T (w) = 



exp 
1 



( huj \ 

\k B T) 



coth 



2k B T 



1 



(144) 



(145) 



Negative frequencies u correspond to the possibility of 
emitting photons into the line 



N T (-\u\) = -N T (\u\)-l. 



(146) 



The Bose-Einstein expression Nt (w) is expected from 
the Hamiltonian of the line, which reads, with the a op- 
erators, 



H = 



E 



+■00 



\cu\ a a [to] a a [—a;] dcj. 



(147) 



We can now give the expression for the anticommutator 
of the fields 



{{aF ± [u)i},a 7± [u 2 ]}) T = 2M T 



U>1 - UJ 2 



S (UJI + LU 2 ) 



wi - uj 2 \ (h{u-i- uj 2 )\ 
sgn [ ) coth [ — — — — — ]d(wi+ lo 2 ) . 



ik B T 



(148) 



(1481 



Equation ( 30 ) with no external drive is identical to Eq. 



, , r i sgn (w) , / huj 



sgn I wj 



N T (H) 



(149) 
(150) 



We now introduce the forward-propagating and 
backward-propagating voltage and current amplitudes 
obeying 

(151) 
(152) 



V~* (x,t) = *JZ C A-* (x,t), 
V*~ (x, *) - v^c^"" («, *) , 



(a:,t) = (x,t)/Z c 
I+- (x,t)=V^ {x,t)/Z c 



(153) 
(154) 



Quantum-mechanically, the voltage and current ampli- 
tudes become hermitian operators 



(x,t) -» (x,t) , 
7^ (x,t) -» t** (x,t) . 



(155) 
(156) 



These operators, in turn, can be expressed in terms of 
field annihilation operators as 



V* (x,t) 



r* (x,t) 



hZc 

47T 



4ttZ c 



+00 



dcu 



+00 



uj\a we 



-iw(f x/vp) 

(157) 
(158) 



All physical operators can be deduced from these primary 
expressions. For instance, the transmission line charge 
operator, describing the charge in the line brought from 
one end to the position x, is 



P (x, t) = i 



+ OO 



do; 



a a; e 



-iuj(t =p x/v p ) 



(159) 



Nyquist model of resistance: semi-infinite 
transmission line 



We now are in a position to deal with the semi-infinite 
line extending from x — to x = 00, whose terminals 
models a resistance R = Z c [see Fig 
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at x 

In that half-line, the left- and right-moving propagat- 
ing waves are no longer independent. We will now re- 
fer to the wave amplitude (x = 0, t) as A in (i) and 
A~* (x = 0, t) as A out (t). The quantum-mechanical volt- 
age across the terminal of the resistance and the current 
flowing into it satisfy the operator relations 



V (t) = V out (t) 

i (t) = i out (t) - 



- V 1D (t) 

i in (t)- 



(160) 
(161) 



These relations can be seen either as continuity equations 
at the interface between the damped circuit and the re- 
sistance/line, or as boundary conditions linking the semi- 
infinite line quantum fields A ln (t) and A out (t). From the 
transmission line relations, 



we obtain 



^ouUn ^ = ^jout.in ^ 



I(t) = Jt V(t)-2I™(t), 



(162) 

(163) 
(164) 
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For a dissipationless circuit with Hamiltonian 
Hbare Qj , where $ is the generalized flux of 
the node electrically connected to the transmission line, 
and Q its canonically conjugate operator (top panel of 
Fig. 16), we can write the Langevin equation, 



dt^ h 

i 

= n 



^barei Q 
H bare j Q 



Rdt 



3> + -= A m (t) . (165) 



The latter equation is just a particular case of the more 
general quantum Langevin equation giving the time evo- 
lution of any operator Y of a system with Hamiltonian 
Hbare, which is coupled to the semi-infinite transmission 
line by an Hamiltonian term proportional to another sys- 
tem operator X , 



dt h 



bare 



Y 



1 

2ik 



X,Y 



, 2R^ 2 A in (t) 



dt 



(166) 



The value of £ in Eq. ( 166[ ) depends on whether the 
damping is "parallel" (( — — 1) or "series" type (C = +1) 
[see Fig. [l6]. In the parallel case, the greater the line 
impedance the smaller the damping, whereas in the series 
case the situation is reversed. 



Equation ( 166 ) should be supplemented by 



and 



i in (h),A in (t 2 ) 



A out (i) = C 



ih d 
T3(*i-t 2 ) 



5{h-t 2 ) (167) 



A in (t)-R c/2 —X 
v ; dt 



(168) 



It follows from the last three equations that the output 
fields have the same commutation relation as the input 
fields 



A out (ti) ,A out (i a ) 



ih 



d 



2 d{h-t 2 ) 



6(h-t 2 ). (169) 



Quantum Langevin equation in the RWA 
approximation 



a narrow range around the oscillator frequency uj a . We 



start from Eq. ( 165 ) and use 

<£ = $ ZPF (a + at) 
(a - at) 



Q = Q 



ZPF 



I 

where $ ZPF = ^/hZj2 and Q ZPF = v / h/2Z a . 

We then obtain, neglecting the effect of driving terms 
oscillating at twice the resonance frequency, 



(170) 
(171) 



df 
with 



0> — "T [Hbare , ft] 



Z a 

" a 2l? a 
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-ni Ain{t) (172) 



A iD (t) 



A iD [uj}e- lult duj. 



(173) 



The field amplitude A ln (t) is non-hermitian and contains 
only the negative frequency component of A ln (t). For 
signals in a narrow band of frequencies around the reso- 
nance frequency, we can make the substitution 



fku, 



A in (t)^a in (t). 



where 



a in (i) 



(174) 



(175) 



The input field operator a m [u;] is identical to a^[w] of 
the infinite line. We finally arrive at the RWA quantum 
Langevin equation, also referred to in the quantum op- 
tics literature as the quantum Langevin equation in the 
Markov approximation 



dt 



[H b . 



where 



7a 



a in (i) , a in (i') f = 6 (t - t') . 



a + y/^ad in {t) 



(176) 



(177) 



For any oscillator, the input output relationship is ob- 
tained from 



/^a(t) = ~a in (t)-Ca out (t) . 



(178) 



We now consider an approximate form of the input- 
output formalism which is valid when the system degree 
of freedom consists of an oscillator with very low damp- 
ing, and for which all the frequencies of interest will lie in 



It is worth noting that although a m and a° nt play the role 
of and in Eq. ( 140 1 , only the average values of 
the moments of a m can be imposed, a out being a "slave" 
of the dynamics of a m , as processed by the oscillator. 
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